Rebuilding cytoskeleton roads: Active-transport-induced polarization of cells 
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Many cellular processes require a polarization axis which generally initially emerges as an inho- 
mogeneous distribution of molecular markers in the cell. We present a simple analytical model of a 
general mechanism of cell polarization taking into account the positive feedback due to the coupled 
dynamics of molecular markers and cytoskeleton filaments. We find that the geometry of the orga- 
nization of cytoskeleton filaments, nucleated on the membrane (e.g. cortical actin) or from a center 
in the cytoplasm (e.g. microtubule asters), dictates whether the system is capable of spontaneous 
polarization or polarizes only in response to external asymmetric signals. Our model also captures 
the main features of recent experiments of cell polarization in two considerably different biological 
systems; namely mating budding yeast and neuron growth cones. 

PACS numbers: 87.10.-e, 83.80.Lz, 87.17.Jj 



Cell polarization is an essential step for many biological 
processes such as cell migration, division, development 
and morphogenesis in widely varying cell types. Polar- 
ization is characterized in its early stages by an inhomo- 
geneous distribution of specific molecular markers. Often 
cell polarization is driven by an external asymmetric sig- 
nal as in the well known example of cells migrating along 
a chemical gradient performing chemotaxis [1] . In mating 
yeast the external signal is a pheromone gradient, which 
causes the cell to grow an elongation known as a shmoo 
in the direction of the pheromone source [1] . In the case 
of nerve cells, external gradients direct the growth of the 
growth cone [2]. However observations show that some 
systems can also polarize spontaneously in the absence of 
any external gradient, e.g. as reported in neutrophil cell 
migration, algae cell division, or mating yeast [3, 4, 5, 6]. 
These two distinct polarization processes, driven or spon- 
taneous, are necessary for cells to fulfill different biolog- 
ical functions, but what determines whether a cell can 
polarize spontaneously or only in response to an external 
signal is still not well understood. In this Rapid Com- 
munication we propose a general physical mechanism to 
address this fundamental biological question, briefly dis- 
cussing budding mating yeast and neuron growth cones 
as specific examples. 

The molecular basis of spontaneous and driven cell po- 
larization has been much debated over the past decade, 
and is likely to involve several pathways including his- 
torical marks (such as that left by the previous budding 
event in yeast [1]). However, it is widely recognized that 
the cytoskeleton plays a crucial role in cell polarization. 
The efficiency of formation of polar caps in yeast is re- 
duced when actin transport is disrupted and the polar 
caps formed are unstable [6, 7, 8]. In the case of neurons, 
it has been shown that the polarization of the growth 
cone is suppressed when microtubules are depolymerized 
[2]. To account for these observations, it is generally ar- 
gued that the cytoskeleton filaments mediate an effective 



positive feedback in the dynamics of polarization markers 
[5, 6, 9]. This arises from the molecular markers not only 
diffusing in the cell cytoplasm, but also being actively 
transported by molecular motors along cytoskeleton fila- 
ments, the dynamic organization of which is regulated by 
the markers themselves (see [10] for another example). 

A number of recent studies have embarked upon mod- 
eling aspects of polarization. Many are reaction-diffusion 
systems in which polarization emerges as a type of Tur- 
ing instability [11, 12, 13, 14, 15, 16] or due to stochastic 
fluctuations [4] and some [5, 6, 16] include cytoskeleton 
proteins as a regulatory factor. In particular Wcdlich- 
Soldncr et al. [6] consider a one-dimensional (ID) model 
which shows that a positive feedback favors the emer- 
gence of polarized states, however, as in other existing 
models, the dynamics of markers in the cytoplasm and 
the geometry of the filaments are not considered. Here we 
propose an analytical model which explicitly addresses 
the geometry of the cytoskeleton and its dynamical cou- 
pling to the transport of markers while remaining gen- 
eral enough to be applicable to a wide range of cell sys- 
tems. Our model provides a minimal mechanism of cell 
polarization induced by active transport and shows that 
the dynamical organization of cytoskeletal filaments ac- 
tually plays a crucial role, since it dictates the polariza- 
tion ability of cells - spontaneous or driven. Specifically, 
our analysis suggests that in general, the case of micro- 
tubule mediated transport does not lead to spontaneous 
polarization but the response to external signals can be 
accurate; while the case of cortical actin mediated trans- 
port leads to spontaneous polarization but a less robust 
response to gradients. 

We model polarization markers as particles which can 
be either on the membrane or in the cytoplasm of the cell. 
For simplicity we assume that the cell is essentially bi- 
dimensional and we neglect curvature effects. The mem- 
brane boundary is then taken as a ID line along the x- 
axis and the cytoplasm is parametrized by (x,z). This 
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geometry is sketched in Fig. 1. The dynamical equa- 
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FIG. 1: Model geometry. Filaments nucleated at the (a) cell 
center and (b) membrane. Arrows indicate the velocity of 
directed transport. Concentration of molecular markers u on 
the membrane and c in the bulk. 

tions for the concentration of markers on the membrane, 
p(x, t), and in the cytoplasm, c(x, z, t), are given by; 

d t p(x, t) =D m d%p(x, t) + j m (1) 
dtc{x,z,t) =- V jb (2) 

where j m ,& are the flux of molecules arriving at the mem- 
brane and the flux in the bulk respectively; 



j m =k on c(x,0,t) - k oS p(x,t) 
jb = D b Vc(x, z, t) + vc(x, z, t). 



(3) 
(4) 



D m ,b are the diffusion constants in the membrane and 
the bulk of the cell respectively and fc on ' off are the rates 
of attachment and detachment to the membrane. In the 
bulk we assume that the rate of molecules switching be- 
tween filament bound and free is high such that the flux 
of particles in the cell bulk, jb, is effectively directed dif- 
fusion where v is the velocity of directed transport of 
markers along cytoskeleton filaments [17, 18]. 

The directed transport v, which drives the system out 
of equilibrium, depends on the geometry of filaments, 
which we assume to be determined by the concentration 
of markers on the membrane ij,(x, t) . This makes the term 
involving v in Eq. 2 a nonlinear coupling term. We con- 
sider two types of idealized geometries for active trans- 
port motivated by the two general biological classes of (a) 
microtubule cytoskeleton and (b) cortical actin systems. 

(a) Filaments that grow from a nucleating center in 
the cytoplasm towards the membrane can be taken in a 
first approximation to be perpendicular to the membrane 
as sketched in Fig. la. Assuming a "search and capture 
model" [19], their orientation is regulated by the polar- 
ity markers due to preferential stabilization and growth 
in regions of high marker concentration. With this ge- 
ometry, in our model, the velocity of directed transport 
can be written simply as v = — ap,[x, t)u z where a is a 
parameter for the strength of coupling. In general mi- 
crotubules arc found with such geometry, e.g. in neuron 



growth cones, for which observations suggest that the po- 
larity markers, GABA receptors, are associated with and 
regulate the growing microtubule ends [2, 20]. 

(b) Filaments that are polymerized from nucleators 
localized at the membrane can be approximated by a 
superposition of asters centered on the membrane as 
sketched in Fig. lb. We model this as a field at posi- 
tion r = (x, z) proportional to the filament concentration 
which decreases with distance from the nucleation point 
r' = (a/, z') on the membrane such that: 

fL/2 „ _ 



r'|2 



p(x', t)dx' 



(5) 



-l/2 \r- 

where a is the coupling parameter and L is the cell 
perimeter. Note that we expect that other decaying func- 
tions of r — r'| would not qualitatively change the re- 
sults. This geometry mimics the organization of corti- 
cal actin (see also [21]), for instance in budding yeast 
[6, 22] in which molecular markers (e.g. Cdc42, Spa2, 
septins) are transported along the filaments towards the 
membrane by myosin V molecular motors. When active, 
these molecules in turn induce actin nucleation at the 
membrane thereby creating a positive feedback. 

Finally, the above equations are completed by the con- 
dition for the conservation of total number of polarity 
markers: M = J^ L 2 . 2 dx p(x, t) + dx J" dz c(x, z, t). 
Note that we can assume that the system is not bounded 
in the z direction since we will only consider concentra- 
tion profiles exponentially decaying with z. 

Spontaneous polarization. To determine whether 
spontaneous polarization occurs in such systems we per- 
form a linear stability analysis of Eqs. (1-2) around 
the homogeneous out of equilibrium steady state solu- 
tion which reads fi° = |^c°(0), c°(z) = c°(0)e _Az . Here 
A = ap°/D b for case (a) and A = apPn/Db for case 
(b). The conservation of markers fixes the value of c°(0) 
in terms M. We perturb the system in the infinite size 
L limit about the homogeneous steady state such that 
fi{x,t) = p° + ti k e ikx+st , c(x,z,t) = c Q (z) + c k {z)e lkx+s \ 
and obtain to linear order from (1): 

fc on c fc (0) 



/'A- 



s + D m k 2 + k ofC ' 



(0) 



In case (a) of a schematic microtubule system where 
v = — afi{x, t)u z , Eq. (2) is solved by Ck(z) = Ae~ Xz + 
(c k (0)-A)e-P z with p = |±|(A 2 + 4fc 2 + 4s/£» 6 ) 1/2 and 
where A is determined by the conservation of flux at the 
membrane (j m + jb ' u 2 )| 2= o = 0. Substituting this into 
(2) leads to the dispersion relation for s(k); 



(s + D b k 2 ) ( fc off (2A - p) + (s + D m k 2 ){\ ~P--^) 



+ X^D b k oS (X -p)=0. 



(T) 



In the limits of both small and large k this gives the same 
solution s = —Dk 2 where D = min(Z?h, D m ). These lim- 
iting negative solutions suggest that the real part of s(k) 
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is always negative, which can indeed be checked by a nu- 
merical search of the parameter space. This shows that 
the homogeneous state is linearly stable and there is no 
spontaneous polarization. This actually compares well to 
experiments in neuron growth cones [2]. When the gra- 
dient was removed, the distribution of receptors returned 
to a symmetric distribution indicating the stability of the 
homogeneous state and the lack of spontaneous polariza- 
tion, as expected from our results. 
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FIG. 2: Solution s(fc) of Eq. (8). k and A are plotted in di- 
mensionless units normalized by the inverse of the cell radius, 
= 2-k/L. Parameter values used are: k° n = lp,ms~ , 
k° s = 0.1s -1 , D b = O.l/WV 1 , D m = 0.01^m 2 s _1 and 
7? = 10/im. The dotted line shows the small k asymptote 



We now consider case (b) of a schematic cortical actin 
system modeled by Eq. (5). In this case Eq. (2) is solved 
by cfe(z) = Be-( x+k *> z + (c fe (0) - B)e~P z , where B is 
determined by the conservation of flux at the membrane 
and the dispersion relation finally reads: 

/ l.on 

(s - \D b k) k oS (2X -p) + (s + D m k 2 )(X - 



D b 



+ X 2 D b k oS (X + k - p) = 0. 



(8) 



When the coupling to the cytoskeleton is switched on 
(A > 0), we find that to linear order in k, Eq. (8) gives 
the positive solution s = D j^ a \k\, which indicates that 
the homogeneous solution is unstable and spontaneous 
polarization occurs. The full solution s(k) is plotted in 
Fig. 2, revealing a maximum at finite k = fc max which cor- 
responds to an instability of finite characteristic length 
27r/fc max . Fig. 3a shows that the most unstable mode 
kmax increases with A, and therefore with the strength 
of the coupling a and the number of markers M. In 
the case of a cell of perimeter L, we conclude that for a 
weak coupling 27r/fc max > L and no patches will be seen. 
Increasing the coupling a (or M), A reaches a series of 
thresholds defined by 27r/fc max = L/n p and above which 
n p patches on the cell will grow. This generic behavior 
predicted in the case of cortical actin mediated trans- 
port agrees qualitatively with available experimental ob- 
servations of spontaneous polarization in yeast. Indeed, 




FIG. 3: a) The most unstable mode A; max as a function of 
A = afj, n/Df, (both normalized by R 1 as in Fig. 2). Car- 
toons of the distribution of markers are shown above the 
curve, b) r = l/s(fc max ) as a function A. The calculated 
solid line is compared to data (points) from [9] for the mea- 
sured polarization time t against pheromone concentration, 
c p h- We used a scaling of A = 2c p h and r = 0.14i to obtain 
the observed fit. Other parameter values used are as in Fig. 2. 



it is found in [9] that spontaneous polarization occurs 
only above a threshold value of uniform concentration of 
pheromone (InM a-factor, which is assumed to drive the 
coupling a of our model), while it is shown in [6] that 
the number of patches indeed grows with the quantity of 
markers (Cdc42). Note that our model also predicts that 
increasing the size of the cell would increase the number 
of patches, which could be tested experimentally. Finally, 
the variation in the timescale for the polarization of a cell, 
t = l/s, with the parameter A is shown in Fig. 3b and 
corresponds qualitatively well to the data from [9] show- 
ing decreasing time for spontaneous yeast polarization 
with increasing concentrations of pheromone. 

Driven polarization. We now show that our model 
also accounts quantitatively for driven polarization, for 
instance in the presence of an activator gradient. We 
assume that such a gradient results in an asymmetric 
activation of the markers on the membrane and there- 
fore in a spatial inhomogencity in the coupling parameter 
a. We consider a perturbation a = oq + ake lkx (where 
k = 2tt / L mimics a cell in a constant gradient) and there- 
fore look for solutions at linear order /i(a;) = /i° + Hke lkx , 
c(x,z) =c°{z) + c k (z)e lkx . 

In case (a) of a microtubule system, Eq. (2) with this 
perturbation is solved by Ck{z) = Ade~ XdZ + (cfc(O) — 
A d )e-P« z with p d = % + i(A^ + 4fc 2 ) 1 / 2 , where X d = 
aopP/Di, and A d is fixed by conservation of flux at the 
membrane. This solution gives the polarized out of equi- 
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librium steady state of the cell in response to an external 
gradient (see Fig. 4), and qualitatively reproduces ex- 
periments on neuron growth cones [2] . Note that such a 
response is linear in the perturbation: it therefore occurs 
for arbitrary small activator gradients and no threshold 
is involved. Furthermore, the response is linearly stable, 
which implies that fluctuations in the external activator 
concentration will be damped out in the response, and 
the cell will polarize on average only along the activator 
gradient. 




FIG. 4: Density of polarity markers c(x, z) for driven po- 
larization with a(x) (a) plotted in the {x,z} model plane 
with parameters as in Fig. 2 and A = 6, k = 1, M = 1000, 
ctk/cto = 0.2 (b) mapped onto a cartoon circular cell. Darker 
shading represents higher density. 

With the same method wc can find the stationary state 
also in case (b) of a cortical actin system. This analy- 
sis however only applies in the regime of small coupling 
ao were the steady state solution is stable as discussed 
earlier. In the unstable regime in an activator gradient 
the cell is likely to polarize spontaneously in the 'wrong' 
direction due to concentration fluctuations, with a prob- 
ability that increases for a short polarization time r (sec 
Fig. 3) . There is therefore a trade off for the biological sys- 
tem between the abilities to polarize fast and to robustly 
polarize in the direction of a weak gradient. Sensitivity 
is compromised in cases such as mating yeast where it is 
important to polarize fast to gain a mating partner. We 
expect this to be the case for many short range cell-cell 
interactions where local gradients are strong, masking 
the effects of fluctuations thereby decreasing the chance 
of spontaneous polarization in the wrong direction. On 
the other hand in cases such as neuron growth cones, the 
speed of polarization is sacrificed in favor of an accurate 
response to weak gradients. 

To conclude, we have shown that the positive feedback 
triggered by the coupled dynamics of molecular polarity 
markers and the cytoskeleton is not sufficient to produce 
spontaneous polarization. We find that whether the sys- 
tem can polarize spontaneously or not actually depends 
crucially on the geometry of cytoskeletal filaments. Ac- 
tive transport of markers along filaments oriented in a 
centered aster, such as microtubules, leads to a stable 
homogeneous state, and polarized states occur only in 



response to a spatial gradient of activator concentration. 
This result is compatible with [4], where it was argued 
on the basis of a generic one dimensional model with lo- 
cal positive feedback only that homogeneous states are 
indeed stable, and that polarized states occur only tran- 
siently due to stochastic fluctuations. On the contrary, 
in the case of filament asters centered on the membrane, 
such as cortical actin, we find that active transport in- 
duces a non local coupling which destabilizes the homo- 
geneous states and therefore leads to spontaneous polar- 
ization above a threshold uniform concentration of acti- 
vator. Since the robustness of polarization along a weak 
gradient is compromised if the cell is able to polarize 
spontaneously this is of great importance in biological 
situations in which either robust response to a gradient 
or fast spontaneous polarization is (dis)advantageous. 
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